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A. 


ABSTRACT 

A  fiber-reinforced  ceramic  subject  to  tensile  stress  in  the  fiber 
direction  can  undergo  extensive  matrix  cracking  normal  to  the  fibers,  while  the 
fibers  remain  intact.  In  this  paper,  the  critical  conditions  for  the  onset  of 
widespread  matrix  cracking  are  studied  analytically  on  the  basis  of  fracture- 
mechanics  theory.  Two  distinct  situations  concerning  the  fiber-matrix 
interface  are  contemplated:  (i)  unbonded  fibers  initially  held  in  the  matrix 
by  thermal  or  other  strain  mismatches,  but  susceptible  to  frictional  slip, 
and  (il)  fibers  that  initially  are  weakly  bonded  to  the  matrix,  but  may  be 
debonded  by  the  stresses  near  the  tip  of  an  advancing  matrix  crack.  The 
results  generalize  those  of  the  Aveston-Cooper-Kelley  theory  for  case  (1) . 
Optimal  thermal  strain  mismatches  for  maximum  cracking  strength  are  studied. 


and  theoretical  results  are  compared  with  experimental  data  for  a  SiC  fiber, 
llthium-alumlna-slllcate  glass  matrix  composite, 
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INTRODUCTION 

Fiber-reinforced  ceramic  materials  have  promising  potential  for  high- 
temperature  applications  (Prewo  and  Brennan,  1980).  Under  tensile  loading  of 
the  composite  in  the  fiber  direction,  the  brittle  matrix  can  undergo  extensive 
cracking  normal  to  the  fibers,  but  the  associated  matrix  cracking  stress  may  be 
substantially  greater  than  the  catastrophic  fracture  stress  of  the  unreinforced 
ceramic.  Furthermore,  with  the  fibers  intact,  the  composite  material  can 
continue  to  sustain  additional  load  up  to  the  fiber-bundle  fracture  stress. 

This  behavior  is  illustrated  by  the  schematic  stress-strain  curve  shown  in 
Fig.  1.  The  slope  of  the  initial  straight  portion  of  the  curve  is  closely 
approximated  by  the  rule  of  mixtures  based  on  matrix  and  fiber  moduli.  Extensive 
matrl)  cracking,  often  involving  a  small  stress  drop,  occurs  at  A  ,  and  the 
matrix  becomes  permeated  by  many,  more-or-less  equally-spaced  cracks  that  traverse 
the  full  cross-section  of  the  specimen.  Under  continued  loading,  the  fibers 
alone  provide  most  of  the  subsequent  stiffness.  The  ultimate  strength  would 
Ideally  be  associated  with  fracture  of  uniformly  strong  fibers,  but  in  practice 
is  degraded  somewhat  as  fibers  fracture  sequentially  rather  than  simultaneously 


before  the  peak  stress  at  B  is  reached. 

In  this  paper,  critical  conditions  for  the  onset  of  widespread  matrix 
cracking  are  studied  theoretically  on  the  basis  of  fracture-mechanics  theory. 

Two  distinct  situations  concerning  the  fiber-matrix  interfaces  are  considered: 

(i)  unbonded  fibers  held  in  the  matrix  by  Initial  pressures  due  to  thermal  or 
other  strain  mismatches,  but  susceptible  to  frictional. slip,  and  (il)  fibers 
under  initial  radial  tension  that  are  weakly  bonded  to  the  matrix,  but  may  be 
debonded  by  the  high  stresses  near  the  tip  of  an  advancing  matrix  crack. 

The  study  of  case  (i) ,  which  generalizes  the  well-known  Aveston-Cooper-Kelly 
(ACK)  theory,  (Aveston  et  al.,  1971  ;  Aveston  and  Kelly,  1973;  Kelly,  1976; 
Aveston  and  Kelly,  1980;  Haxmant  et  al.,  1983)  is  based  on  the  analysis  of 
steady-state  crack  growth  in  the  matrix.  The  concept  adopted  (slightly  different 
from  that  of  ACK)  is  that  a  "first"  planar  crack  will  propagate  across  the 
composite  under  an  applied  stress  that  becomes  constant  during  the  propagation 
as  soon  as  the  crack  engulfs  more  than  a  few  fibers.  With  dynamic  effects 
neglected,  the  stress  associated  with  this  steady-state  cracking  is  equivalent 
to  the  "first  cracking"  stress  of  ACK.  (The  initiation  of  growth  of  the  most 
critical  flaw  in  the  matrix  could  require  a  somewhat  higher  stress  than  that 
associated  with  steady-state  growth  —  hence  the  slight  dip  in  the  stress-strain 
curve  of  Fig.  1.)  Figure  2  illustrates  the  matrix  crack  as  it  proceeds  across 
the  composite.  With  enough  frictional  resistance,  no  slip  will  occur  at 
the  interfaces,  as  shown  in  Fig.  2(a).  When  slip  does  occur  (Fig.  2(b)),  the 
slip  length  along  the  fibers  on  either  side  of  the  crack  can  be  expected  to 
approach  an  asymptotic  value  on  the  downstream  side  of  the  crack  front. 

The  presumption  that  Coulomb  friction  provides  the  resistance  to  fiber 
slip  implies  that  positive  fiber-matrix  pressures  are  imposed  by  strain 


mismatches  that  occur  during  the  fabrication  process.  It  does  not  however 
follow  that  increasing  such  mismatches  would  necessarily  raise  the  matrix  cracking 
stress,  despite  the  larger  frictional  resistance  thereby  provided.  The  same 
strcin  mismatches  also  generally  lead  to  Initial  axial  tensile  stresses  which  act 
to  reduce  the  cracking  strength.  Accordingly,  optimal  strain  mismatches  can  be 
expected  to  exist,  and  these  will  be  estimated. 

In  the  case  of  bonded  fibers,  the  matrix  cracking  strength  will  depend  on 
the  debonding  toughness  of  the  interface.  In  the  presence  of  sufficiently  high 
debonding  toughness,  the  first  matrix  crack  will  propagate  in  a  manner 
indistinguishable  from  that  of  the  no-slip  frictional  case  (Fig.  2(a)).  If 
debonding  does  occur,  and  the  Interface  pressure  is  negative  (i.e.,  tensile 
residual  stresses  exist  between  fiber  and  matrix)  the  debonded  regions  will  open 
up,  and  the  crack  will  propagate  as  shown  in  Fig.  2(c).  The  steady-state 
cracking  calculation  will  be  made  for  this  case  on  the  basis  of  an  elementary 
analysis  of  the  debonding  process  near  the  advancing  crack  front. 

ENERGY  RELATIONS 

A  fairly  general  relation  will  be  derived  for  the  loss  in  potential  energy 
of  a  prestressed  elastic  body,  within  which,  under  constant  additional  load, 
cracks  develop  and  open  up,  and  also  sliding  occurs  along  internal  interfaces. 
These  relatione  will  be  used  subsequently  in  the  steady-state  cracking 
calculations.  Figure  3  shows  three  states  of  the  body.  In  state  (0),  the  body 
is  free  of  external  load,  but  contains  an  initial  tensor  stress  distribution 
00  in  its  volume  V  .  With  the  external  vector  tractions  T  applied  to  the 
external  boundary  in  state  (1) ,  the  stress  becomes  ,  and  additional 

displacements  u^  ,  compatible  with  additional  strains  are  produced.  The 
body  may  now  contain  open  cracks,  as  well  as  Internal  surfaces  in  which  sliding 


has  occurred.  In  state  (2) ,  with  no  change  in  T  ,  more  open  cracking  has 
occurred,  and  additional  frictional  sliding  has  taken  place  along  the  Interface 
Sp  .  The  final  stresses  are  now  a2  *  aQd  the  displacements  and  strains,  still 
measured  from  the  ground  state  (0) ,  are  and  •  We  want  to  calculate 
the  potential  energy  loss  associated  with  the  transition  from  the 

energy  in  state  (1)  to  in  state  (2) .  The  elastic  constitutive 
relations  may  be  written 

e  -  M (o-cr0)  (1] 


for  the  strain  changes  produced  by  o  ,  where  M  is  a  linear  operator.  Then 
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In  each  case,  the  volume  integral  represents  the  elastic  energy  stored  in  the 
body.  Since  C?1  :  M(a2)  -  o2  :  M(ox)  ,  the  energy  loss  may  be  written  as 


iTj-  *2  -  \  :  M(°l"°2)dV  "  |  T*<Vu2)dS  (5' 

V  Sj 

We  now  assume  that  in  state  (2),  the  shear  tractions  on  S?  act  in  a  direction 
opposite  to  that  of  the  relative  sliding,  and  have  constant  magnitude  Tg  . 
Then,  by  the  principle  of  virtual  work, 

[T*(u,-u„)dS  -  fo„  :  M(0,-0„)dV  -  T  [  |Av|dS 


(6 


where  |Av|  is  the  magnitude  of  the  relative  slip  on  Sj,  that  has  occurred 
during  the  transition  from  state  (1)  to  state  (2).  Hence 

"l"1^  "  2  J(ai“°2)  :  M<°1-°2)dV  +  t8  J  lAv'dS  (7) 

v  sF 

If  we  assume  further  that  slip  on  has  been  unvarying  and  monotonic  in 
direction  during  the  transition  to  state  (2) ,  and  that  the  sliding  resistance 
has  always  been  equal  to  T  ,  the  frictional  energy  £_  dissipated  (as  heat) 
is  precisely  the  last  integral  in  (7).  Hence 

*2  "  \  }(crr°2)  :  <ere2)dV  +  5F  (8) 

V 

under  the  stipulated  assumptions.  This  result  is  clearly  not  valid  under 
conditions  of  variable-direction  slip,  or  history-dependent  frictional 
resistance,  during  the  transition  from  state  (1)  to  state  (2).  It  does  remain 

correct  for  polntwise  variations  in  t  . 

8 

STEADY-STATE  CRACKING  RELATIONS 

To  apply  the  energy  relations  just  developed,  we  contemplate  a  long  matrix 
crack  of  length  s  in  a  very  wide  specimen  of  width  W  and  unit  thickness,  as 
shown  in  Fig.  4.  The  crack  extends  through  the  thickness  of  the  specimen,  with 
a  straight  front  CC  ,  but  all  of  the  fibers  are  intact.  With  no  change  in  the 
average  applied  stress  o  ,  the  crack  is  presumed  to  advance  an  amount  As  to 
C'C’  .  Now  identify  the  initial  uncracked  and  unloaded  but  possibly  prestressed 
state  of  the  specimen  with  state  (0) ,  and  let  the  states  before  and  after  the 
crack  advance  As  correspond  to  states  (1)  and  (2).  The  assumption  of  steady- 
state  cracking  means  that  the  stresses  at  the  crack  front,  averaged  through  the 
thickness,  remain  unchanged  during  the  crack  growth,  and  also  that  the  upstream  and 


downstream  stress  states  far  ahead  of  and  behind  the  crack  front  do  not 
change.  Consequently,  if  we  define  Py  and  as  the  upstream  and  downstream 
potential  energies  per  unit  cross  sectional  area  of  the  composite,  it  follows 
that 

fl_7r2  “  (PlfVA8  (9) 

Hence,  the  potential  energy  release  rate  (per  unit  crack  extension,  per  unit 
thickness)  is 

L  H 

V  PD  "  2A  |  {  !  ^EtTeD^ dAdz  +  ds 

C-LA 

c 

where  ,  Ey  and  oQ  ,  are  the  upstream  and  downstream  stress  and  strain 
distributions,  and  A  is  a  representative  cross-sectional  area  of  the  composite. 

H.  c 

Here  -gjj-  is  the  frictional  energy  dissipation  rate  (per  unit  thickness) 
associated  with  fiber-matrix  slip. 

In  the  case  of  unbonded,  frictionally  constrained.  Blipping  fibers,  the 

energy-release  rate  Py  -  Py  must  be  balanced  by  the  sum  of  this  frictional 

energy  dissipation  rate  and  the  critical  matrix  crack-extension  energy-release 

rate  c  t  per  unit  thickness  of  the  composite.  Hence  the  relation 
m  m 

L 

■dr  J  I  <W  ■  (11) 

C  -L  A 

c 

governs  matrix  cracking,  for  both  the  slip  and  no-slip  cases.  If  slip  does 
occur,  the  validity  of  this  result  requires  that  there  be  monotonlcally 
increasing  slip  along  each  fiber. 

In  the  case  of  initially  bonded,  debonding  fibers,  the  frictional  term  in 
(10)  is  absent,  but  now  a  debonding  energy  release  absorbs  part  of  P.  -  P_  . 


-8- 


For  a  unit  crack  advance,  the  increment  in  debonded  surface  area  per  fiber  is 

2ita on  each  side  of  the  crack,  and  the  number  per  unit  area  of  newly  debonded 
2 

fibers  is  c^/(na  )  ,  so  that  the  debonding  energy  release  rate  is  4cf (i^/a) , 
Hence 

L 

2A~  \  \  (0U"°D)  :  (eu'£D)dAdz  *  Cm*m  +  4cf<Va>*d  (12) 

C  -L  A 

c 

In  order  to  implement  (11)  and  (12)  for  the  calculation  of  the  cracking 
stress,  we  now  have  to  estimate  0^  and  0^  ,  and  for  the  debonding  situation 
we  also  must  estimate  the  debonding  length  2^  . 

FIBER-MATRIX  STRESS  ANALYSIS 


Upstream  Stresses 

Far  ahead  of  the  crack  tip,  the  axial  stresses  in  the  fibers  and  the  matrix 
in  the  loaded  composite  are  those  of  the  uncracked  material.  The  upstream 
stresses  are  therefore  well  approximated  by 


o"  -  (Ef/E)0  +  0* 


(13) 


oJJ  -  (E  /E)a  +  o*  ] 
in  in  m  / 

where  a  is  the  average  applied  stress,  o*  and  0*  are  the  initial  axial 

stresses  in  the  unloaded  composite,  and  Ee  ,  E  are  the  fiber  and  matrix 

t  m 

Young's  moduli.  This  approximation  neglects  the  effects  of  transverse  stresses 

on  axial  strains  ,  and  is  consistent  with  the  rule-of-mixtures  expression 

E  -  c.E,  +  c  E  (14) 

ft  mm 

for  the  effective  axial  modulus  of  the  composite.  The  initial  and  total  stresses 
satisfy 


C.oI  +  C  0* 

f  f  mm 


0 


(15) 


-9- 


respectlvely. 


c,0  +  C  o  ■  o 

f  f  mm 


Downstream  Stresses 


Behind  the  crack  tip,  the  average  fiber  and  matrix  axial  stresses  at  the 


crack  face  are 


°£  m  0/ct 


) 

and  for  L»a  the  stresses  at  z«L  are  given  by  Eqs.  (13).  Approximate  , 
shear-lag  analyses  will  provide  the  far-downstream  stress  distributions  in 
each  of  the  cases  shown  In  Fig.  2. 

No-slip  case 

Far  from  the  crack  tip,  an  Isolated  composite-cylinder  shear-lag  model, 
similar  to  that  adopted  by  Aveston  and  Kelly  (1973)  will  be  used.  Each  fiber  is 
presumed  to  be  embedded  in  a  matrix  cylinder  of  outer  radius  R  chosen  as 

R  -  a/^  (18) 

to  provide  the  correct  volume  concentration  of  fibers  (Fig.  5(a)).  The  model  is 
further  simplified  (Fig.  5(b))  by  concentrating  all  of  the  axial-stress-carrying 
area  of  the  matrix  at  an  effective  radius  R  between  a  and  R  ,  and  assuming 
that  the  region  in  a<r<R  supports  only  shear  stresses  T  (r,z)  .  The 
equilibrium  and  constitutive  relations  in  this  region  simplify  to 

3t  t 

-5r  +  -r-°  (19> 

V  ■  G.  If  <20> 

where  w(r,z)  is  the  axial  displacement,  measured  from  the  uncracked  state. 

It  follows  that 


where  T^(z)  is  the  interface  shear  stress,  given  by 


,  v  Gm(VWf} 

x  (Z)  -  - (22) 

a  log(R/a) 

in  terms  of  the  fiber  and  matrix  displacements  w.  ■w(a,z)  and  w  *w(R,z)  . 

£  in 

Fiber  equilibrium  implies 

3  0.  - 

-5r+y,i‘°  <23> 

and  since  the  composite  cylinder  is  isolated,  Eq.  (16)  remains  valid. 


Eliminating  w^  and  w^  from 


Of -Of  dwf  ^ 
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0-0  dw 
m  _ m  m 


(24) 


E  dz  ' 

m 

and  Eqs.  (16),  (22),  (23),  and  applying  the  boundary  conditions  (17)  and  (13) 
at  z»0  and  j  z  |  ■« «»  respectively  leads  to  the  following  results  for  the 
downstream  stresses: 

D  U  .  .  U  -plzj/a  ^ 
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where 


2G  E 
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If  the  fibers  are  held  in  the  matrix  by  friction,  this  result  is  valid  only  if 


-li¬ 


the  no-sllp  condition 


% a  TJ<0+) 
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is  met. 

Slipping  fibers 

When  the  no-slip  condition  (27)  is  violated,  frictional  sliding  between  fiber 

and  matrix,  with  T.  ■  T  ,  occurs  in  a  length  £  on  either  side  of  the  crack. 

X  8  8 

Then  Eqs.  (16),  (17)  and  (23)  imply  that 

o°  -  o/cf  -  2Tg|z|/a 

<£-2<cf'c m)Tslzi/a  ►  (28) 

) 

i  8 

for  0<S  |z|^£g  ,  far  downstream  from  the  crack  tip.  Re-solving  Eqs.  (16), 
(21)-(23)  ,  now  with  the  boundary  conditions  T.  ■  T  at  |z|-£  and,  again, 

X  8  8 

(13)  at  |  z  |  *  00  ,  leads  to 
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for  |z|i£g  .  Then  imposition  of  the  requirement  that  the  axial  stresses  be 
continuous  at  |z|*£  provides  the  equation 


a+  (E/E 
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APPENDIX  A 


Estimate  of  Effective  Radius  R  in  Shear-Lag  Model 
Evaluation  of  the  shearing  energy  contribution  in  Eq.  (32)  gives 


2c,(l+v  ) 
i  in 

E 
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log  R/a 


00 

I  Tid2 
•00 


(Al) 


This  relates  to  the  simplified  model  of  Fig.  5(b),  but  now  let  us  return  to  the 
configuration  of  Fig.  5(a),  and  contemplate  a  continuous  shear  lag  model  in  which 
the  matrix  stress  0^  is  distributed  across  the  outer  cylinder.  Longitudinal 
equilibrium  Implies 


3a  3(rr  ) 
_ m  ,  _ rz 

r  3z  3r 


(A2) 


together  with  the  boundary  conditions 


Trz(a)  *  Ti 
Tra<R)  -  0 


(A3) 


The  assumed  distribution 


rz  a  VCJ  r 


(A4) 


satisfies  these  boundary  conditions.  Substitution  into  (A2)  gives 


3o 

-gf  -  2(Cf/cm)(Tl/.)  (A5) 

uniform  in  r  ,  and  this  is  consistent  with  the  equilibrium  requirements  (23)  and 
(16).  Accordingly,  an  appropriate  solution  for  the  z-distributions  of  the  stresses 
based  on  the  principle  of  minimum  complementary  energy  (which  requires  the  use  of 
an  equilibrium  approximation)  based  on  the  assumption  (A4)  is  legitimate.  Except 
for  the  definition  of  the  characteristic  parameter  p  ,  this  solution  will  be 
identical  to  the  one  derived  on  the  basis  of  the  simplified  model.  The 


propagation  of  matrix  cracks  into  the  fibers.  The  post-matrix-cracking  strength 
exhibited  in  Fig.  1  and  the  accompanying  pseudo-ductility,  would  then  be  lost. 
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Measured  values  of  a cr  gave 

(®cr)exp  "  290420  ^ 

These  results  suggest  the  presence  of  a  small  initial  axial  compression. 
This  is  consistent  with  observations  reported  by  Marshall  and  Evans  (1985)  of 
matrix  crack  closure,  upon  unloading,  at  small  tensile  loads.  But  initial  axial 
compression  in  the  matrix  would  ordinarily  be  accompanied  by  tension  normal  to 
the  fiber-matrix  interface.  Accordingly,  interfacial  roughness  rather  than 
Coulomb  friction  may  have  been  the  primary  source  of  the  Interface  shear 
resistance  in  the  Marshall-Evans  experiments. 

CONCLUDING  REMARKS 

The  two  idealized  assumptions  pursued  herein  concerning  the  fiber-matrix 
Interface  —  frlctlonally  constrained,  sliding  fibers  and  initially  bonded, 
debonding  fibers  —  are  not,  of  course,  exhaustive.  Combinations  of  these 
possibilities  could  coexist,  and  Interface  roughness  might  play  a  more  important 
role  than  interface  pressure  in  providing  slipping  resistance.  Two  interesting, 
if  tentative,  conclusions  can  nevertheless  be  reached: 

1.  If  Coulomb  friction  is  operative,  optimal  strain  mismatches  exist 
that  maximize  the  matrix  cracking  strength. 

2.  In  the  case  of  initially  bonded  fibers,  a  fairly  small  interface 

debonding  toughness  suffices  to  inhibit  debonding 

during  the  matrix  cracking  process. 

A  final  cautionary  note:  the  inhibition  of  either  debonding  or  slipping 
may  be  quite  undesirable  despite  the  fact  that  the  matrix  cracking  strength  is 
thereby  Increased.  Rill  maintenance  of  fiber-matrix  continuity  facilitates 


and  this  condition  will  generally  be  easily  met.  This  means  that  as  a  function 

of  decreasing  (4^/J^)  ,  the  matrix  cracking  strength  will  drop  abruptly  from 

its  no-debond  value  when  {.t  ,1b  )  falls  below  its  debond  threshold  value.  This 

d  m 

is  illustrated  in  Fig.  12  for  Er/E  -3  ,  V  -1/4  ,  and  o1-©  .  It  is 

t  m  m  m 

interesting  to  note  that  once  it  drops  from  its  no-debond  value,  O  remains 

cr 

fairly  insensitive  to  debond  toughness  until  extremely  low  values  of  b,  are 

d 

reached . 


EXPERIMENTS 


In  some  recent  experiments,  Marshall  and  Evans  (1985)  studied  first-cracking 
in  a  ceramic  system  consisting  of  silicon-carbide  fibers  in  a  lithium-alumino¬ 
silicate  glass  matrix.  The  nominal  values  of  pertinent  parameters  were 
cf  -  .5 
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Push-through  and  Indentation  tests  of  individual  fibers  in  composite  samples  gave 

measured  values  of  T  *2.0  MPa  ,  suggesting  that  the  frictional-slip  model  should 

8 

be  applicable. 

On  the  basis  of  the  nominal  data,  Eqs.  ( 38) — (40)  give 

aQ  -  1625  MPa  (B-.88) 
a1  -  265  MPa 


This  put 8  a1f°Qm  *16  well  within  the  large-slip  range  (Fig.  7),  and  so  the 
theoretical  prediction  for  the  cracking  stress  is 
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axis  and  each  lower  branch  must  be  overcome.  Since  the  debonded  lengths  along 
the  lower  branches  are  quite  small,  it  seems  reasonable  to  presume  the  presence 
of  initial  material  flaws  and  imperfections  that  are  equivalent  to  initially 
debonded  regions  of  similarly  small  size.  Then  crack-tip  stresses  would  push 
these  effectively  debonded  lengths  to  the  upper  branches  of  the  curves. 

Threshold  Bond  Toughness 

The  function  Q(X)  in  Fig.  9  attains  a  maximum  value  Q* »  2.061  at 

*  * 

X  -  .9204.  Since  debonding  can  not  occur  for  Q>Q  ,  this  critical  value  of 

Q  can  be  used  in  Eq.  (62b)  to  obtain,  as  a  function  of  c^  ,  threshold  values 

0t  Jb  )*  of  the  ratio  of  debonding  to  fracture  toughness  that  would  prevent 
d  m 

debonding.  The  results,  shown  in  Fig.  11,  indicate  that  a  debonding  toughness 
that  is  quite  small  in  comparison  to  matrix  fracture  toughness  suffices  to 
suppress  crack-tip  debonding  over  the  practical  range  of  fiber  fraction. 

Critical  Cracking  Stress 


With  the  use  of  Eq.  (36)  for  p  ,  Eq.  (55)  for  Ocr  becomes 


!sE  +  xi. 

°o  E*  °0 


B2r  6e  i1/2f*di 

cmL(1+Vm)EfJ  J 


where  B  (Fig.  6)  is  given  by  Eq.  (39)  as  a  function  of  c^  ,  and  £^/a  is 

defined  as  a  function  of  c.  and  JSJb  by  Eqs.  (62)  and  the  upper  branch  of 

t  am 

the  curve  in  Fig.  10.  For  UtJb  )  >  ibjb)  ,  no  debonding  occurs,  and  the 
s  a  m  a  m 


a  a 

old  no-slip  result  ~  r  +  g~  — ■  1  applies.  When  debonding  does  occur,  we  can 

0  m  0 

expect  the  result  (64)  to  provide  a  lower  cracking  strength  than  that  for  the 
no-debond  case.  For  <  U&lbJ*  »  the  right-hand  side  of  (64)  is  certainly 


less  than  unity  if 


°vV<  £[ 


f>-  Efa+v. 
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The  solution  o£  (59) ,  with  the  boundary  conditions  U( ±  i .)  -  0  ,  permits  the 

d 

calculation  of  the  energy  expression 


(z)U(z)dz 


which  represents  the  loss  in  potential  energy  of  half  of  the  loaded  matrix  jacket 
due  to  debonding.  Then  the  energy-release-rate  relation 

3V 

**bt  ‘  TT  <61> 

a 

provides  the  condition  governing  l ^  .  The  result  (Appendix  C)  is  the  pair  of 


parametric  equations 


a  -  (1  + 
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relating  l.l a  to  3} .lb  ,  where 
a  am 


is  plotted  in  Fig.  9. 
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Note  that  Eqs.  (62)  imply  that  S^/a  is  not  a  single-valued  function  of 

b Jh  .  But  from  the  energy-based  derivation  that  led  to  this  result,  it  follows 
d  m 

that  combinations  of  A./a  and  H  Jit  associated  with  the  region  to  the  left 

a  dm 

of  the  curve  in  Fig.  9  are  unstable,  while  those  to  the  right  are  stable.  The 
implications  of  this  are  better  seen  in  Fig.  10,  which  show,  as  examples, 
explicit  plots  of  i^/a  vs.  for  several  values  of  c^  .  The  upper 

part  of  each  curve  represents  stable  debonded  lengths,  but  for  debonding  to 
occur  at  all  an  energy  barrier  associated  with  the  region  between  the  horizontal 


•*.  •*.  **«  *V**«  **  •  ^ 


We  suppose  that  stress  changes  in  the  vicinity  of  the  crack  tip  debond 

the  fiber-matrix  interface  just  ahead  of  the  advancing  crack  for  a  distance  l. 

d 

on  each  side  of  the  plane  of  the  crack.  To  estimate  we  will  again  adopt 
a  composite  cylinder  model  in  which  we  pretend  that  axisymmetric  debonding  is 
produced  by  an  axisymmetric  distribution  of  load  applied  to  the  matrix  cylinder 
that  jackets  the  fiber  (Fig.  8).  Conservatively,  we  take  the  magnitude  of  this 
loading  as  the  stress 

K 

0  (z)  .  - 2_  (5( 

and  apply  it  at  the  mean  radius  R  of  the  matrix  jacket.  Here 

ID 

2  1/2 

K  -  [E  H  / (1— v  )]  is  the  critical  elastic  Btress- intensity  factor  of  the 
n  q  n 

matrix,  and  (56)  is  just  the  asymptotic  distribution  of  horizontal  tension  just 
above  the  crack  tip.  The  debond  length  ftg  will  be  calculated  on  the  basis  of 
an  energy  balance  involving  the  debonding  toughness  % ^  ,  and  the  energy 
changes  in  the  matrix  during  debonding.  The  deformation  of  the  matrix  will  be 
analyzed  on  the  basis  of  thick-cylinder  theory,  in  which  only  transverse  shear 
stresses  and  circumferential  tension  resist  radial  displacement.  In  terms  of 
the  cylinder  thickness 

-1/2 

t  ■  R  -  a  -  a(c^  -1)  (5‘ 

and  the  mean  radius 

\ -  i<R+a)  ■  mi/2+i>  (s 

the  differential  equation  governing  the  radial  displacement  U(z)  at  r“R 


is  taken  as 


linear  initial  part  of  the  curve  in  Fig.  7,  justifying  the  use  of  Eq.  (40)  for 


There  is  an  interesting  connection  between  and  the  mismatch 

for  self-cracking  of  the  matrix  in  the  absence  of  external  stress.  Under  the 
assumption  of  large  slip,  self-cracking  occurs  when  the  right-hand  side  of  (47) 
vanishes,  which  Implies  that 

fiSC  “  3^  R0PT 

MATRIX  CRACKING:  (li)  INITIALLY  BONDED,  DEBONDING  FIBERS 


Cracking  Condition 

We  assume  that  there  is  initial  transverse  Interface  tension  in  the  composite, 
so  that  downstream  from  the  crack  tip,  regions  of  the  fiber-matrix  interface  that 
have  been  debonded  stay  open. 

The  steady-state  cracking  condition  (12)  is  then 
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(54) 


where  we  will  use  (13),  (31)  for  |z|  <  ,  and  (25),  with  |z|  replaced  by 

(|z|-l^)  ,  for  |z|  >  .  (As  before,  we  ignore  transverse-stress  contributions 

to  strain  energy.)  This  leads  to  the  result 


where 

length 
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fd  fd 
a  it. 
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is  given  by  (36),  but  now  we  require  an  estimate  for  the  debond 


(55) 


are  obtained. 


If  we  anticipate  that  the  optimum  value  of  ft  will  lead  to  large-slip 
matrix  cracking,  it  is  appropriate  to  use  Eq.  (40)  for  0cr  ,  with  Ts*yq 


in  Eq.  (41).  Then 


a  [6yc*E  *11/3f  'i  1/3  o' 
cr  i  f  m  I  q  |  m 
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with  q/E  and  o  / E  given  by  (45).  For  ft>0  ,  0  _  will  attain  a  maximum 

m  mm  cr 

value  at  ft*  ^0PT  when  the  condition 


is  met .  Hence 
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At  this  optimum  design,  the  associated  slip  length,  based  on  the  first  term  of 


Eq.  (30),  is  given  by 


3X,E. 

V*  ’  nsr 


Corroboration  of  the  validity  of  the  large-slip  assumption  follows  from  the 


observation  that 


<°„  w  (E/v^ .  ir.“)1/2FE(ltyf 

»k  J  L  6K-  j 


at  ft-ft^^.,  .  This  value  will  generally  be  small  enough  to  fall  in  the  nearly 


Is  valid  for  the  sliding  frictional  resistance  Tg  ,  it  follows  that  increasing 

the  strain  mismatch  would  raise  0.  as  well  as  0*  .  Then  Eq.  (40)  implies  the 

i  m 

existence  of  an  optimal  strain  that  maximizes  the  externally  applied  cracking 

stress  0  .  (However,  the  assumption  of  Coulomb  friction  is  not  necessarily 

valid.  Conceivably,  Tg  may  be  due  primarily  to  interface  roughness,  in  which 

case  the  optimization  study  that  follows  is  inapplicable.) 

Suppose  that  non-elastic  strains  e  and  e,  occur  isotropically  in  matrix 

in  i 

and  fiber  during  fabrication,  and  call 

n  -  (ef-ea) 

the  strain  mismatch.  (If,  for  example,  the  mismatch  is  due  only  to  thermal 

strains,  Q*  (a. -a  )AT  ,  where  AT  is  the  temperature  change,  and  a  ,  a,  are 
cm  mi 

the  linear  thermal  expansions  over  the  range  AT  .  Note  that  0  is  positive 
if  the  matrix  contracts  more  than  the  fiber  during  fabrication.) 

For  simplicity,  we  assume  that  both  the  fibers  and  the  matrix  are  isotropic. 
A  straightforward  analysis  of  the  composite  cylinder  model  of  Fig.  5(a)  then 
gives 


where  X,  ,  X,  are  functions  of  c.  ,  E,/E_  ,  v.  ,  and  v  shown  explicitly 
12  z  I  m  t  m 

in  Appendix  B.  For  cf-*-l  ,  X^X^l  ,  and  both  ^  and  X2  will  not  vary 
too  much  from  unity  for  reasonable  values  of  c^  and  E^/E^  *  ^  we  assume 


v  ■  v  -  V  ,  the  neat  forms 
f  m 
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X  -  1  -  -kl-E/E,)  -  4(1+E/E,) 


(46) 


the  results  can  be  manipulated  Into  the  form 


where 


(43) 


In  the  range  Y>1  specified  by  the  slip  condition  (27),  Bqs.  (43)  are  parametric 
relations  giving 


(o  +E/E  )0 
cr  m  m 


as  the  function  of  the  independent  variable  (O^/Oq)  plotted  in  Fig.  7.  For 
Y«1  ,  (01/aQ)»31/3- 1.442  ,  and  so  for  (o^Oq)  >  31/3  ,  the  no-slip  result 


o  +  (E/E  >0*" 
cr _ m  m 


-  1 


applies.  For  Y-*-«*  ,  (o^Oq)  0  ,  and  °cr+ (E/Em)o^  approaches  the  large-slip 
ACK  value  .  This  large-slip  result  is  a  good  approximation  over  the 
substantial  initial  portion  of  the  curve  in  Fig.  7  that  is  nearly  linear.  In 
the  slipping  range,  the  slip  length  lg  may  be  found  from  Eq.  (30). 

Optimal  Strain  Mismatch 

Mismatches  between  the  non-elastic  fiber  and  matrix  strains  that  occur 
during  fabrication  (e.g.,  due  to  cooling,  plasticity,  creep,  or  phase  trans¬ 
formation)  will  produce  Initial  matrix  stresses  o*  as  well  as  fiber-matrix 
interface  pressure  q  ,  and  a  positive  q  will  generally  go  along  with  positive 
o1  .  If  a  Coulomb  friction  lav  of  the  form 


Then  B-*l  for  c^l  >  and,  as  shown  in  Fig.  6,  B  does  not  vary  much  from 

unity  over  a  large  range  of  fiber  concentration  c^  . 

The  ACK  results  for  slipping  fibers  (Aveston  et  al.,  1972)  can  be  recovered 

from  Eq.  (32)  by  substituting  Eqs.  (28)  for  0?  and  oD  in  |z| <  £  ;  dropping 

r  m  s 

the  shear  contribution  to  the  energy;  neglecting  all  energy  contributions  in 

|z|  >  l  ;  and  dropping  the  (1/p)  term  from  the  formula  (30)  for  l  / a  .  (It 
8  8 

can  be  verified  that  these  truncations  are  all  asymptotically  valid  for  1&/ 

The  result  is 


(40) 


where 


is  equivalent  to  the  ACK  expression  for  the  critical,  large-slip  cracking  strain. 

(Aveston  et  al.  do  not  actually  present  a  counterpart  to  Eq.  (40);  they  derive 

the  cracking  criterion  0  /E»Oj/E  for  o*-0  ,  and  separately,  also  deduce 

that  a^/E-o, /E  is  the  condition  for  self-cracking  in  the  absence  of  external 
m  i 

loading.) 

To  bridge  between  the  no-slip  and  large-  lip  ACK  results,  we  can  substitute 

the  full  expressions  (29)  for  |  z  |  ^  £  ,  together  with  (28)  for  |z|*Jl  . 

s  s 


into  (32).  With  Y  defined  by 
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wherein,  by  Eq.  (21),  T®  - -t®  .  (Here  we  neglect  contributions  to  strain  energy 

rz  it  1 

associated  with  downstream  changes  in  transverse  stress.)  For  the  no-slip  case 
substitution  of  Eqs.  (25)  into  (32)  leads  to  the  formula 


0  a1  [c.E.XfpT 

E  EB  L  aEBEJ 


for  the  cracking  stress  a  .  Except  for  the  initial-stress  term  cj/E  , 

cr  m  m 

this  is  essentially  the  result  originally  given  by  Aveston  and  Kelly  (1973)  for 
no-slip  matrix  cracking. 

Aveston  and  Kelly  do  not  specify  R  in  the  definition  (26)  for  p  ,  beyond 
the  unelaborated  statement  that  R  is  equal  to  the  radius  at  which  the  matrix 
displacement  equals  its  average  value.  On  a  different  basis,  the  explicit 
estimate 

2  log  c.+  c  (3-c,) 

log  R/a  -  - - 5 -  (34) 

in 

is  derived  in  Appendix  A.  (This  gives  (R-a)/(R-a)  *►  1/3  for  1  ;  for 
-  -3/4 

Cj-*0  ,  (R-a)/(R-a)  e  *  .47  .)  If  we  introduce  the  utility  constant 


6  log  R/a 


into  the  definition  of  p  ,  we  get 


B2f  6E  V 

c,Uf«+vJ 


and  we  can  rewrite  (33)  in  the  convenient  form 
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for  the  slip  length.  For  £.g“0  this  Is  consistent  with  the  no-slip  requirement 

(27),  and  Eqs.  (28),  (29)  are,  of  course,  valid  only  for  l  /atO  .  For  It 

8  8 

larger  than  a  few  fiber  radii,  the  contribution  to  fcg/ a  of  the  term  (1/p)  in 
(30)  will  be  small. 

Initially  bonded,  debonding  fibers 

Now  we  suppose  that  debonding  along  a  length  l.  on  either  side  of  the  crack 

a 

is  produced  near  the  crack  tip  by  interface  fracture,  and  that  the  debonded 
regions  remain  open  downstream  of  the  crack  tip.  The  axial  and  shear  stresses  in 
0£  |z|  £  are  simply 


0D  ■  0 
m 


(31) 


and  for  |z|£Jtd  shear-lag  analysis  reproduces  the  results  (25)  from  the  no-slip 
case,  with  |z|  in  the  exponents  simply  replaced  by  (|z|-  K^)  .  Ue  presume 
that  the  interface  shear  stresses  will  not  produce  any  additional  debonding  that 
increases  1,  beyond  its  crack-tip  value. 

U 

Note  that  in  contrast  to  the  case  of  frictionally  slipping  fibers,  the  shear- 
lag  solution  for  debonded  fibers  involves  a  discontinuity  in  at  |z|»i^  . 

MATRIX  CRACKING:  (i)  UNBONDED,  FRICTIONALLY  CONSTRAINED  FIBERS 

Critical  Cracking  Stress 

Using  the  stresses  of  the  shear-lag  model  in  the  steady-state  cracking 


equation  (11)  gives  (for  L -*■<») 


correspondence  between  the  two  solutions  is  easily  found  by  replacing  R  by 


R  in  the  shear  energy  contribution  in  Eq.  (32)  and  using  (A4)  for  T 
This  gives 


rz 


+v  )  fCr\  rf  ^2  *  ^2  .  e  \  f  <i 

m  I  f|  |R|  ,  R,,  3|R]  lfal  1  f  2 

~k)  Ira  log;+1-*ra  -*fei  J  J  v 


2(l+v_)  (C 
E 


2 

With  (R/a)  -  (l/cf)  ,  comparison  of  (Al)  and  (A6)  then  provides  the  estimate 
(34)  for  R/a  . 


APPENDIX  B 


Initial  Stresses 

By  the  classical  Lane  solution  for  the  unloaded  composite  cylinder  of 
Fig.  5(a),  the  initial  circumferential  stresses  Og^  and  Og^  in  the  fiber 
and  matrix  at  their  interface  are 


°6m  “  q(1+Cf)/cm  j 

where  q  is  the  interface  pressure.  The  conditions  of  interface  strain 
continuity  are 

+  •»  -  r,0r  V'ef-rt  +  ef  I 

m  f 

ftae»,-VV<l)1  +  *.  -  ri0ef-V0r’>i +  v 

in  x 

Substitution  of  (Bl) ,  elimination  of  O*  via  Eq.  (15),  and  solution  for  o* 
and  q  gives  the  results  (45)  of  the  text,  with  X^  and  Xj  defined  by 

1-  (l-E/Ef)(l-Vf)/2  +  co(Va-vf)/2-  (E/Ef)[vf  +  (Vm-Vf)cfEf/E]2 

X1  (I^TU  “ 

m 

(1-  (l-E/Ef)/2] (1+Vf)  +  (l+cf) (Vjj-Vj)/2 

X2  "  I 

where 

A  ■  1  +  vf  +  (Vn-vf)cfEf/E 


and  E  is  given  by  (14).  Setting  v -  v.  gives  Eqs.  (46). 


APPENDIX  C 


Debond-Length  Analysis 


Letting 


z  -  4d(C/X) 


■-jpefSter. 


(Cl) 


where 


,1/2 


x  ■  <V*>  (137) 

reduces  the  differential  equation  (59)  to 

-u?c  +  u  -  ur1/2 

and  the  boundary  conditions  become  u(  ±  X)  ■  0  .  The  energy-release-rate 
condition  (61)  becomes 

0 

For  £>0  ,  the  solution  of  Eq.  (C3)  is 

C  5 

u.  f  8inh(C-C')dCt  .  cosh  C  f  sinh(X-E')dC' 

J  CO.h  X  J  =r 


(C2) 


(C3) 


(C4) 


(C5) 


which  gives 


5u  m  cosh  C 
3X  "  cosh2X  J  JV 


A 

f  CQ8h  C'dC* 
J  J7T 


(C6) 


Then  (C4)  provides  the  results  (62),  (63)  of  the  text,  when  c^  is  Introduced 
via  Eqs.  (57),  (58). 

For  calculation  purposes,  the  function  Q(X)  in  Eq.  (63)  may  be  written  in 


terms  of  Dawson's  Integral 


D(Z)  ■  e  J  e  ds 
0 

which  is  tabulated  and  available  in  software.  Thus 


Q(X)  - 


!XD(X1/2)+-5  erf(X1/2) 


fc)+-y-  erf < 
cosh  X 


Z  2 

where  erf(Z)w  —  /  e  8  ds  .  For  X  large,  Q~l/X 


(c)  Initially  bonded,  debonding  fibers 


store  iu) 


state  (1) 


stote  (2) 


Fig.  3  Three  successive  states:  (0)  unloaded,  prestressed,  (1)  loaded 
(2)  new  open  crack  surface  S  ,  more  sliding  on  surface  S_  . 


t  — H  K~ 


versus  debond- toughness /matrix' 

o  (v  -  1/4)  . 


no  debonding 


debonding 


.2  .4  .6 

Cf 

.8 

Fig.  11  Threshold  bond  toughness 

(v  -1/4)  . 

m 

^  no  debond 

Fig.  12  Matrix  cracking  atrangtta  veraua  debond  toughness/matrix  toughness 
ratio;  illustrative  example  E./E  •  3  ,  v  ■ 1/4  ,  o*«0  . 


END 

FILMED 

7-85 


DTIC 


